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Abstract. The cellular automaton model for traffic flow exhibits a jam- 
ming transition from a free-flow phase to a congested phase. In the de- 
terministic case this transition corresponds to a critical point with di- 
verging correlation length. We present data from numerical simulations 
which suggest the absence of critical behavior in the presence of noise. 
The transition of the deterministic case is smeared out and one only 
observes the remnants of the critical point. 



1 Introduction 

The cellular automaton (CA) approach to traffic flow theory [Q has attracted 
much interest in recent years. Compared to the earlier attempts in modeling 
traffic flow (see e.g. |Q and references therein) CA models can be used very 
efficiently for computer simulations. This allows to perform realtime simulations 
of large networks |^ . 

One of the simplest models of that type is the model introduced by Nagel and 
Schreckenberg (NaSch model) 0|. This model is discrete in space and time. A 
discrete space means that the road is divided into L cells which can be occupied 
by a car or can be empty. 

The model contains only a few parameters: The braking noise p, the maxi- 
mum velocity v^ax and, if one considers periodic boundary conditions, the aver- 
age density p of cars. The update is divided into four steps which are applied in 
parallel to all cars. The first step (Rl) is an acceleration step. The velocities vj of 
each car j = 1, . . . , iV not already propagating with the maximum velocity Vmax 
are increased by one. The second step (R2) is designed to avoid accidents. If a car 
has dj empty cells in front of it and its velocity (after step (Rl)) exceeds dj the 
velocity is reduced to dj. Up to now the dynamics is completely deterministic. 
Noise is introduced via the randomisation step (R3). Here the velocities of mov- 
ing cars {vj > 1) are decreased by one with probability p. The steps (R1)-(R3) 
give the new velocity Vj for each car j. In the last step of the update procedure 
the positions Xj of the cars are shifted by Vj cells (R4) to xj + Vj . 

For a given set of parameters Vmax and p one can distinguish two differ- 
ent phases depending on the density For low densities each car can 
propagate with its desired velocity Vmax — P such that the fiow is given by 
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J = p{vmax — p)- For higher values of the average density jams occur and start 
and stop waves dominate the dynamics. Here we study the nature of the transi- 
tion and the behaviour of the model at the transition point. 

2 Order - Parameter 

For a proper description of the transition one has to introduce an order-parameter 
(see also with a qualitative different behaviour in the two phases. Due to 
particle conservation one cannot use the density of cars, which would be the ana- 
logue to the magnetisation in the Ising model. Among several possible choices 
of an order-parameter the density of nearest neighbour pairs 

1 ^ 

m{p) = X! niUi+i > (1) 
1=1 

with Tii — 1(0) if the site i is occupied (empty), is probably the simplest choice 
of a local order-parameter. One could think of alternative definitions, e.g. the 
density of cars with a distance less than Vmax , but those results are in qualitative 
agreement with results obtained using (|^). Fig. [l] shows the density dependence 
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Fig. 1. Behavior of the order-parameter for a finite braking probability p. It 
does not vanish exactly for p < pc, but converges smoothly to zero even for 
small values of the braking probability p. 

of the order-parameter for Vmax — 2 and p — 0.50. In the free flow regime the 
order-parameter vanishes and one observes a monotonous increase of m for higher 
values of the average density. At the transition regime m converges smoothly to 
zero if a finite value of the braking noise is considered. This behaviour is not the 
result of a finite system size, the data shown in Fig. |] show no size dependence 
anymore. Due to the smooth convergence of the order-parameter it is difficult 
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to determine a transition density, but it is obvious that interactions between the 
cars occur aheady below the density of maximum flow. 

In the deterministic hmit {p = 0) the situation is different. Here one observes 



a sharp transition at pc = — 



.+1 



the density of maximum flow. 



3 Spatial Correlations 

A striking feature of a second-order phase transition is the diverging length 
scale at criticality. In order to determine the length scale for a given density, we 
measured the density-density correlation function 



G(r) 
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Fig. H shows the behaviour of the correlation function near the transition density 




Fig. 2. Correlation function in the presence of noise {vmax = 2,p = j^g). The 
amplitude of the correlation function decays exponentially for all values of p. 



for Vmax — 2 and p — . The correlation function is a damped oscillation with 
period Vmax + 1- The amplitude is decaying exponentially for all values of the 
average density if a finite value of p is considered. In the deterministic limit one 
gets precisely at p = pc the following behaviour of the correlation function: 



G(r) = 



pI for r = mod(u,„Q2; -|- 1) 
else 



(3) 



Due to the exponential decay one can determine the density dependence of the 
correlation length ^(p) if p > 0. 

Fig. H shows the behaviour of £,{p) for Vmax = 2 and different values of p. The 
figure shows that the correlation length has a maximum ^max near the transition 
density, which is less pronounced for higher values of the braking noise. 

In Fig. ^ the noise dependence of ^max{p) is shown. The numerical data are 
in excellent agreement with 



p^O. 



(4) 
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^(p) , v„„ = 2 



p = 1/128 ^ 
p= 1/64 
p= 1/32 
p= 1/16 «^ 
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Fig. 3. Density dependence of the correlation length in the vicinity of the tran- 
sition. 




Fig. 4. Noise dependence of ^max for different maximum velocities. Independent 
of the maximum velocity, ^maxip) ~ l/\/P holds in the limit p — > 0. 



The divergence for p = is consistent with the behaviour of G(r) at the 
transition density, where one observes a constant amplitude of G{r) and therefore 
an infinite correlation length. Finally it should be noted that the numerical 
results can be confirmed analytically for Vmax = 1 • 

4 Summary 

We found evidence for the absence of criticality in the NaSch model in the 
presence of noise. For finite p the second order transition of the deterministic case 
is smeared out, similar to the situation of a equilibrium second order transition 
in a finite system. Here this effect is caused by the presence of noise, p > 0. 

For small values of p one finds an ordering transition close to pc = ^/{vmax + 
1). The ordering is due the parallel dynamics and can be found for all values of 
the maximum velocity, including Vmax = 1- Larger values of the noise p favour the 
formation of jams and a tendency to phase separation occurs (see also [^[ |To| , pl| ) 
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if Vmax > 1- The demixing of free flow and the jammed regime becomes more and 
more exact if larger values of Vmax and p are considered, because interactions 
between cars become rare in the free flow regime of the system. 

The effect of fluctuations is twofold: First, one already finds blocked cars 
for densities below pc and second, exact phase separated states are not stable, 
e.g. one observes the spontaneous formation of jams in the outflow region of a 
megajam. 

Finally we discuss the deterministic limit p ~ 1. Here wc find mctastable 
states for densities below pc = l/3 ii Vmax > 1- These mctastable states consist 
of moving cars with at least two empty cells in front. They are mctastable in 
the following sense: If a car stops due to a pertubation of the system, all other 
cars have to stop after 0{L) lattice updates. These mctastable states can only 
be found for Vmax > 1, because for Vmax = 1 obviously the cars cannot move. 
The lack of mctastable states for < p < 1 is consistent with our observation 
that exact phase separation is absent for this cases. 
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